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A  Gibbs-type  partition  of  the  chaotic  motion  is  introduced.  The  distribution  function  of  the  spec¬ 
tral  density  defined  on  such  partition  is  Gaussian.  The  variance  of  this  distribution  is  the  Fourier 
transform  of  the  correlation  function.  We  demonstrate  this  by  direct  numerical  computations  for 
the  simple  models  of  chaos.  These  results  are  the  consequence  of  translational  invariance  and 
should  be  valid  for  the  general  case  of  chaotic  motion  described  by  differential  equations. 


The  purpose  of  this  paper  is  to  establish  a  general  rela¬ 
tionship  for  fluctuation  of  the  spectral  density  of  the 
chaotic  motion,  which  is  similar  to  the  Einstein  fluctua¬ 
tion  formula  in  statistical  mechanics.1 

We  begin  with  the  simplest  model  of  chaotic  motion 
given  by  one-dimensional  quadratic  mapping2 

x,  + 1  =  ox,(  1  — x, )  .  (1) 


Here,  t  =0,1,2, .  .  .  can  be  thought  as  a  discrete  time. 
The  iterative  sequence  produced  by  Eq.  (1), 
j x0,xl,x2,  ■  ■  .  ,xr ),  will  be  called  an  orbit;  T  is  the  time 
length  of  the  orbit  and  x0  its  initial  condition.  It  is  easy 
to  see  that  for  the  values  of  parameter  0<a<4,  the  or¬ 
bits  are  bounded,  0<x,  <1.  The  chaotic  orbits  are 
characterized  by  a  positive  Liapunov  number  X  >  0:1 


X=  lim 


t  =o 


dx , 


(2) 


On  the  other  hand,  in  the  case  of  regular  (periodic)  orbits 
X<0.  Another  important  property  of  chaotic  orbits  is 
their  power  spectrum.  The  Fourier  transform  for  a 
discrete  variable  x,  can  be  written  in  the  form 


1  -r- 

1  t=  0 

1  V 

X'=V¥  ^a"€ 

1  n  =0 


(3) 

(4) 


Here,  co„  =  2trn  /T,  a_n  —a*.  Consider  the  limit  T  —  oc. 
In  the  case  of  chaotic  orbits,  the  to  spectrum  becomes 
continuous  and  we  can  substitute 


r-i 


rI,m  T  2  - 

T^cc  1  -1=0 


1 

2n 


l 


2ir 


dco  . 


(5) 


In  the  case  of  periodic  orbits,  a(w)  is  discrete  (in  the  limit 
T— *  oo  ),  that  is,  a  sum  of  delta  functions.  Often,  chaotic 
motion  also  exhibits  delta-functional  peaks  in  the  spec¬ 
trum.  For  example,  usually,  a0  =  {x,)V T  — 6(e>!  when 
T —*  oo.  Such  delta-functional  peaks  could  also  happen 
for  a  finite  value  of  to. 3 

Introduce  now  a  correlation  function: 


C(r)=  (x,  +  Tx, )  -  (x,  )2  . 


(6) 


The  time  averaging  is  defined  by  the  equation 


(  A  (x, ))  =  lim 

T  —  ao 


l  r-' 

~  2  AU,) 
1  1  =0 


(7) 


Here  A  (x)  is  some  function  of  x.  The  correlation  func¬ 
tion  C(t,T)  averaged  over  the  finite  interval  T  is  simply 
related  to  power  spectrum.  Using  Eqs.  (4),  (6),  and  (7)  we 
can  get 

Cir,T)=~T±' \alcon)\V^r  .  (8) 

1  n  =  1 

Notice  that  the  term  a0  =  (x,  )VT  is  absent  here  due  to 
the  subtraction  in  Eq.  (6).  We  have  done  detailed  numer¬ 
ical  studies  of  the  correlation  function  C(r,T)  and  the 
power  spectrum  |  a(w„ )  ] 2  for  simple  models  of  chaotic 
motion.  The  typical  behavior  of  C(r,  T)  is  plotted  in  Fig. 
1  for  a  quadratic  mapping.  The  value  of  the  parameter 
a  =  3. 94  corresponds  to  a  chaotic  attractor  with  a  posi¬ 
tive  Liapunov  number  X  =  0.54.  In  Fig.  1,  T —  2 000. 
First  of  all,  there  exists  numerically  a  well-defined  limit 
C (t,T  —<■  cc  )  —  C{t).  Long-time  behavior  of  C(r— «•  x  )  in 
many  cases  can  be  approximated  as  an  oscillation  and  an 
exponential  decay  with  characteristic  time  rc. 

The  plot  of  a  power  spectrum  j  a(con  )  1 2  for  the  same 
value  of  parameters  is  given  in  Fig.  2.  Due  to  the 
discreteness  of  our  model  and  the  symmetry  condition 
a(  —  o)n  )  —  at(a>n  )  we  have  plotted  it  only  in  the  interval 


FIG.  1.  Correlation  function  C(r )  computed  for  quadratic 
mapping,  Eq.  (1).  The  values  of  parameters  are  or  -  *  -M 
r  =  2000,  and  X  =  0.54. 
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0  <  id  <  tt.  It  is  obvious  from  this  plot  that  a  I  con )  is  a  fluc¬ 
tuating  quantity.  A  small  change  in  the  initial  condition 
of  the  orbit  *0  or  *n  the  value  of  two  adjacent  frequencies 
2 vn /T,  2ir(n  +  l)/T  will  usually  result  in  an  abrupt 
change  in  the  value  of  a(con).  Most  importantly,  there 
exists  no  limit  a(a>, ,  T— *  x  ).  Thus,  the  left-hand  side  of 
Eq.  (8)  is  well-defined  in  the  limit  T  — » -x  but  each  term  in 
the  sum  on  the  right-hand  side  is  not.  From  the 
mathematical  point  of  view  we  have  do  some  resumma¬ 
tion  in  the  right-hand  side  of  Eq.  (8)  in  order  to  have  an 
average  quantity  {  !a(«„)j2)  which  converges  in  the 
limit  T— -  oc  to  the  Fourier  transform  of  the  correlation 
function.  From  the  physics  point  of  view  this  means 
choosing  an  ensemble  over  which  such  averaging  can  be 
done.  We  will  introduce  now  such  an  ensemble.  Let  us 
use  as  a  new  time  of  averaging  MT  instead  of  T  in  Eq.  (8); 


here  M  »  1  is  a  very  large  integer.  We  also  assume  that 
T  »tc.  Now  we  split  the  time  averaging  in  Eq.  (7)  into 
two  parts 


<  A  >  w  I  I 


Tm 

2 

■  Tim- 


A(x,  )=<  (  A  )T  >( 


(9) 


Here,  the  brackets  with  subscript  T  denote  usual  time 
averaging,  while  brackets  with  the  subscript  G  denote 
averaging  over  the  ensemble  of  M  test  orbits.  We  call 
this  a  Gibbs-type  ensemble  because  in  the  case  of  a  one¬ 
dimensional  lattice  when  oilr  discrete  time  is  changed  to 
a  discrete  spacing,  it  is  equivalent  to  a  Gibbs  partition  of 
statistical  mechanics.1  Then,  we  can  rewrite  Eq.  (8)  in 
the  form 


T 


C(t,MT)=  -zz  2  (\a(ojn)\2)c 


(10) 


We  have  computed  numerically  the  value  of  {  |  a(w,  II2) 
(the  subscript  G  will  be  omitted  from  now  on)  as  a  func- 
*ion  of  a>„  for  the  example  given  in  the  previous  Figs.  1 
and  2(a).  It  is  plotted  in  Fig.  2(b)  for  a  =  3. 94,  T  =  128, 
and  M  =  105.  The  average  power  spectrum  converges  nu¬ 
merically  to  a  well-defined  continuous  function  of  cu  in 
the  limit  M —*  oo,  T—*co.  This  function  is  a  statistical 
property  of  the  orbit  and  is  independent  of  the  initial  con¬ 
dition  x0.  We  should  mention  here  that  in  general 
several  attractors  could  be  present,  then  depending  on  the 
value  of  initial  condition  x0  the  orbit  will  evolve  on  the 
different  attractors  and  could  have  different  values  of 
{  \  a(co„)\2).  Consider  now  the  limit  Af— »oo,  T-*oo. 
Then  we  can  write  Eq.  (8)  using  Eq.  (5)  in  the  form 


C (r)=  — flwG(co)e‘^ciaj  ,  r>0 
2n  jo 

G  —  (  \a(o))\2)=2CiT)e-,at  . 

r=0 


(11) 


(12) 
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FIG.  2.  (a)  Power  spectrum  |a(<u, )  | 2  computed  for  the  same 
case:  a  =  3.94,  7"=128.  (b)  Average  power  spectrum 

(  a ( I  :  2 )  computed  for  the  same  case  a  =  3.94,  T  =  128,  and 
M  =  10'. 


The  singularity  of  G(tu)  in  the  complex  plane  z=eia 
which  is  closest  to  the  unit  circle  |  z  \  =  1  determines  the 
long-time  behavior  of  C(r—  oo  ).  In  the  case  of  a  pole 
z0  =  exp(  —  1  /tc  +  /'<y0),  |  z0  |  <  1,  we  will  have 

C(t— ►  oo  )~cos(w0r+<p0)e  /Tf  .  (13) 

We  will  turn  now  to  the  numerical  computations  of  the 
distribution  function  P(an  )  of  a  Fourier  coefficient  a, 
with  a  given  frequency  a>n  =2irn  /T.  We  compute  this 
distribution  based  on  Gibbs-type  ensemble  of  M  test  or¬ 
bits  according  to 

Pit  C  )=— - — — 

M  Af,  Af, 

.  |4, 

Here,  AM  is  a  number  of  test  orbits  in  the  interval  Af, 
and  Af„.  We  have  done  extensive  numerical  studies  of 
this  distribution  for  many  different  chaotic  attractors  de¬ 
scribed  by  Eq.  (1)  and  for  other  mappings.4  The  typical 
results  of  these  computations  are  presented  in  Figs.  3  .«nd 
4.  Figure  3  is  a  contour  plot  of  the  distribution  function 
P(an).  These  computation^  were  done  for  a  quadr.u-i 
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FIG.  3.  Contour  plots  of  the  distribution  function  P(a„ ) 
computed  for  quadratic  mapping:  a  =  3.94,  T  =  128,  cj=0.785. 


mapping  Eq.  (9)  with  the  values  of  parameters:  a  =  3. 94, 
7*=  128,  6)=0.785,  M=106.  It  is  clear  from  Fig.  3  that 
P[a„ )  is  axisymmetric.  The  dependence  of  P(  |  an  |  )  has 
been  computed  in  Fig.  4.  The  values  of  the  parameters 
for  this  computation  were  the  same  as  for  Fig.  3,  except 
that  oi=  1.23.  The  results  of  these  computations  can  be 
parametrized  in  the  form 


FIG.  4.  Distribution  function  of  the  spectral  density  comput¬ 
ed  numerically  for  the  case  of  quadratic  mapping  a  =  3. 94, 
T  —  128,  &)=  1.23,  and  M  =  106.  Its  functional  form  is  given  by 
Eq.  (15). 


that  the  real  and  imaginary  part  of  an  have  a  Gaussian 
distribution 

P(a„)  =  P(Z„)P(£n)  ,  (18) 


P(  la. I  >= 


2i  an 


<  !J> 


exp 


1  a«  1 2 

/  1  Dt  b  \  ^  r-rr\ 

[ 

(15)  ^  p 

I  2<s2> 

We  would  like  now  to  show  that  Eq.  (15)  is  a  conse¬ 
quence  of  the  translational  invariance  of  our  distribution. 
Let  us  make  a  transformation  x',=xt  +r,  where  r  is  an  ar¬ 


bitrary  integer, 
transformation 


In  Fourier  space  this  corresponds  to 


=  \a„  |  exp[/(<p„  +o„t)] 


(16) 


Distribution  function  P(a„ )  should  be  invariant  under 
such  transformation.  This  implies  that  it  is  independent 
of  phase  <p„ .  The  Fourier  coefficient  for  the  whole  inter¬ 
val  MT  can  be  expressed  through  our  ensemble  as 


l2-*< 


2) 


(19) 


(20) 


Obviously,  Eqs.  ( 1 8)— (20)  are  equivalent  to  Eq.  (15). 
(  \an  j2)  is  a  function  of  n  and  T;  so  far,  we  have  not 
specified  the  value  of  T.  But  in  the  limit  T»ic  we  have 
always  found  numerically  that  it  becomes  a  function  of 
single  parameter  to  and  independent  of  T. 

Due  to  a  general  nature  of  our  arguments,  they  can  be 
directly  applied  to  a  continuous  time,  namely  chaotic 
motion  described  by  differential  equations. 


a(<on,MT)=-J=-  2  |Oe'Vr.  (17) 

Phases  <p ”  in  the  sum  (17)  are  uniformly  distributed. 
Then,  based  on  the  central-limit  theorem  we  can  claim 
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